XOR-Ising Loops and the Gaussian Free Field 
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We find by simulation that the interfaces in the exclusive-or (XOR) of two independent 2D 
Ising spin configurations at the critical temperature form an ensemble of loops that have the same 
distribution as the contour lines of the Gaussian free field, but with the heights of the contours 
spaced y2 times as far apart as they are for the conformal loop ensemble CLE4 or the double dimer 
model on the square lattice. For domains with boundary, various natural boundary conditions for 
the two Ising models correspond to certain boundary heights for the Gaussian free field. 



XOR-ISING LOOPS AND RELATED MODELS 

The double-Ising model consists of two independent 
Ising spin configurations a and r on the same graph. 
The double-Ising model on planar graphs is closely re- 
lated to a variety of statistical physics models, including 
the Ashkin-Teller model @, §, the 8- vertex model || §, 
the Gaussian model |3j-pl, the double-dimer model, and 
some of which we discuss below. 
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The XOR-Ising configuration £ (illustrated in Figure 
is defined by = a v X t v for each vertex v of the graph 
(if the spin values are ±1), or equivalently, = a v © t v 
(if the spin values are false and true), where © denotes 
exclusive-or (XOR). The XOR-Ising spins £ also called 
the polarization of the double-Ising model ||. The in- 
terfaces between +1 spins and —1 spins in an XOR-Ising 
configuration form loops or paths connecting boundary 
points. We find empirically that when the two Ising mod- 
els are at the critical temperature, the XOR-Ising loops 
are closely related to contour lines of a Gaussian free 
field, but that, in a sense that we will make more precise, 
there are fewer XOR-Ising loops than loops in other loop 
models that are related to the Gaussian free field. 





FIG. 1. (Color online) On the left is the critical Ising model 
on a lozenge-shaped domain of order L = 32 with black- white 
boundary conditions. The interface from bottom to top is in 
the scaling limit SLE3 pl| , |l2^ ]. In the middle is the critical 
Ising model with black boundary conditions. The ensemble of 
loops is in the limit CLE 3 0, On the right is the XOR 
of the left and middle spin configurations. The interface from 
bottom to top appears to be SLE 4 v /j_ 1 yj-i m the scaling 
limit. 



XOR-Ising interfaces are (essentially) the same as the 
doublc-dimer model on the Fisher lattice, as we now ex- 
plain. A dimer configuration on a graph is a pairing of the 
vertices such that every vertex is paired with exactly one 
of its neighbors. Fisher showed that the Ising model on 
any graph has the same partition function as dimer con- 
figurations on a related graph |fuj|| . This relation is not 
combinatorial, but for planar graphs there is a combi- 
natorial correspondence between dimers and Ising spins, 
based on the Ising model's low temperature expansion 
0. When the Ising spins are on the triangular lattice, 
the dimer configurations are on a graph consisting of do- 
decagons and triangles, called the Fisher lattice. 

The double-dimer model is formed by superimposing 
two independent random dimer configurations; the result 
is a collection of loops and doubled edges. If one of the 
dimer configurations has two defects (monomers) on the 
boundary, then the double-dimer configuration also has a 
path connecting the defects. If we ignore the route that a 
double-dimer loop (on the Fisher lattice) takes within the 
triangles, the double-dimer loop traces out a path where 
one Ising configuration has aligned spins while the other 
one does not. Then double-dimer loops on the Fisher 
lattice (ignoring doubled edges and the intratriangle por- 
tions of loops) are just the XOR-Ising interfaces. 

It is instructive to compare XOR-Ising loops with the 
double-dimer model on other lattices. On the square lat- 
tice or hexagonal lattice with suitable boundary condi- 
tions, the scaling limit of a doublc-dimer path connect- 
ing two boundary points is conformally invariant jlij and 
is thought to be described by Schramm-Loewner evolu- 
tion with parameter 4 (SLE 4 ) (see e.g., [^0). The 
scaling limit of the whole ensemble of loops is conjec- 
tured to be the conformal loop ensemble with parame- 
ter 4 (CLE4) |ll|. (The conformal loop ensemble CLE K 
is thought to be the scaling limit of the 0(n) model when 
n = — 2cos(47t/k) JXs| . ) Dimer configurations on the 
square and hexagonal lattices have height functions that 
are known to behave like the Gaussian free field in the 
scaling limit , and CLE4 is known to be the scaling 

limit of contour lines of the Gaussian free field where the 
heights of the contours are separated by a certain spacing 
But the techniques for analyzing the double-dimer 
model fl4[ depend in an essential way upon the lattice 
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quantity 



XOR-Ising loops nearby exact value 



CLE4 loops (exact value) 



loop dimension 1. 5000 ±0.0003 

winding angle variance / log L 1.0000 ± 0.0003 

E[# loops around point] / log L 0.07168 ± 0.00005 

Var[# loops around point] / log L 0.05057 ± 0.00003 

gasket dimension 1.91400 ± 0.00006 
electrical thickness shown in Fig 

difference in log conformal radii shown in Fig 



1.5000- ■• =3/2 

1.0000 - ■■ = 1 

0.07164- ■■ = 1/(V2tt 2 ) 

0.05066- •• = 1/(2tt 2 ) 

1.91421 ••• = 1/2 + V2 
excursion time of BM from [— n, n] 
exit time of BM from \—n, V2tt] 



1.5000- ■■ = 3/2 

1.0000- ■• = 1 

0.10132--- = 1/tt 2 

0.06755--- = 2/(3tt 2 ) 

1.87500 - - - = 15/8 
excursion time of BM from [— it, n] 
exit time of BM from [— 7r, it] 



TABLE I. Comparison of XOR-Ising loops with the conformal loop ensemble CLE4. The (first five) measured values are from 
configurations on L x L tori for L = 2 10 ,2 11 ,2 12 . The properties of individual XOR-Ising loops appear to be the same as for 
CLE4 loops, but the two loop ensembles are different. The measured values of the XOR-Ising loops are consistent with contour 
lines of the Gaussian free field, but with height differences between contours being y/2 times as large as for CLE4. 



being bipartite, which the Fisher lattice is not. Dimers 
on the Fisher lattice have no height function, and the 
monomer correlations for the Fisher lattice |^2|, ^3] be- 
have quite differently from monomer correlations on the 
square and hexagonal lattices Q. Nonetheless, as we 
shall see, individual double-dimer loops on the critical 
Fisher lattice appear to have the same limiting behavior 
as loops of the square- or hexagonal-lattice double-dimer 
models, but there are V2 times fewer loops in the Fisher- 
lattice double-dimer model. 

The polarization of the double- Ising model (the XOR- 
Ising configuration) at the critical temperature has been 
identified with a "spin wave operator" of the Gaussian 
free field || (see also {§-§|), and this identification has 
been used to compute Ising spin correlations |p2| , p3| . 
However, the focus was on correlation functions, and it 
is unclear how to use these works to extract information 
about the geometric properties of XOR-Ising interfaces. 

Boutillier and de Tiliere |], [l(| studied the partition 
function of the double- Ising model on "isoradial" graphs, 
which include the hexagonal and square lattices. They 
showed that when the Ising models are critical, the par- 
tition function is identical to the partition function for 
cycle-rooted spanning forests on the same graph. Span- 
ning trees can be transformed into dimers on a related bi- 
partite graph ^(| , which have a height function whose 
scaling limit is the Gaussian free field. However, there is 
no known algorithm for transforming cycle-rooted span- 
ning forests to pairs of Ising spin configurations, so it is 
not clear how to use their results to relate the Gaussian 
free field to XOR-Ising interfaces. 

Sheffield and Werner gave a characterization of the 
conformal loop ensemble CLE K for 8/3 < n < 4. The 
gasket is the set of points not surrounded by any loop 
of CLE K . Sheffield and Werner showed that the CLE K 
gasket is the same as the set of points not surrounded 
by any loop of the Brownian loop soup with intensity 
c = (3k — 8) (6 — k)/{2k) (the central charge) |27 . Since 
c(k — 3) = i and c(k = 4) = 1, the intersection of two in- 
dependent CLE3 gaskets is a CLE 4 gasket. It follows im- 



mediately that the gasket for the OR (rather than XOR) 
of two independent critical Ising models is the same as 
the CLE4 gasket |28| , which is the scaling limit of contour 
lines of the Gaussian free field. It is not clear, however, 
how to relate the OR-Ising loops to the XOR-Ising loops. 



XOR-ISING LOOPS IN THE BULK 

We generated many XOR-Ising configurations on large 
L x L tori and measure d variou s properties of the XOR- 
Ising loop ensemble (see Table l| ) whose values are known 
for CLE K . These include the length of the longest loop 
(to measure the loop dimension |{29| ) , the winding angle 
variance of the longest loop |30], |3lf, the average number 
of loops surrounding a point |32J, |33j , the variance in the 
number of loops surrounding a point J|3|, and the 
number of points surrounded by no loops (to measure 
the gasket dimension |33|, The first two quantities 

are properties of individual loops, and match the values 
for CLE 4 loops. The latter three quantities are properties 
of the ensemble rather than individual loops, and these 
do not match the values for the CLE 4 ensemble. For 
example, the number of XOR-Ising loops surrounding a 
point is smaller by a factor of v2 than the number of 
CLE4 loops surrounding a disk of radius 1/L. 

In view of the various connections between the double- 
Ising model and the Gaussian free field, and the close 
match between the measured properties of individual 
XOR-Ising loops and CLE4 loops, a natural hypothe- 
sis is that the XOR-Ising loops are distributed accord- 
ing to contour lines of the Gaussian free field, but with 
the height spacing between contours larger than what it 
would be for CLE4 loops. To test this, we measured the 
electrical properties of the XOR-Ising loops. 

Consider the surface of a cylinder of length I 3> 1, 
radius 1, and girth lix made of a resistive medium with 
resistance 27r. Suppose there is a sequence of noncon- 
tractible loops C±, . . . , Ck that wind around the short di- 
rection of the cylinder. Let R-(£i) and denote 
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FIG. 2. The electrical thickness of XOR-Ising loops closely 
matches the Brownian excursion time from [— 7r,7r]. 
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FIG. 3. The difference in log-conformal radii of XOR-Ising 
loops closely matches the Brownian exit time from [— w, V^n]. 



the electrical resistances between loop Li and the left 
and right ends of the cylinder respectively. If the loops 
are CLE4 loops, then — R_{d) 

is distributed according to the exit time of a standard 
Brownian motion from the interval [— 7r, tt] 1 33 , [m), and 
in J34| it was predicted that the electrical thickness 
I — (R-(£i) + R + (£i)) of loops Hi far from either end 
of cylinder is distributed according to the excursion time 
of a standard Brownian motion from the interval [— n, n]. 
Schramm and Sheffield proved that CLE4 is the scaling 
limit of the contour lines of the Gaussian free field with 
a certain spacing A between the contour heights pjj. It 
also follows from their work that when the height spacing 
of the contours is xX, then i?_(£j + i) — i£_(£j) is dis- 
tributed according the exit time of standard Brownian 
motion from the asymmetric interval [— ir,xir] BU p8[. 



with SLE 4 yj-i V2-1 ( see figure 5), which corresponds 
to boundary heights of ±\/2A for the Gaussian free field 
plf . If instead the boundary spins alternate between 
black and white (effectively free boundary conditions), 
with an interface forced by two phase shifts (see Fig- 
ure ||), then the left-crossing probabilities are consistent 
with SLE4,_i _i, which corresponds to boundary heights 
of ±e (in the limit e I 0) for the Gaussian free field. 

We also tried a third set of boundary conditions, which 
specifies the boundary spin values of neither Ising model, 
but instead conditions on their forming an interface. We 
did this by generating Ising spin configurations on a long 
cylinder with antiperiodic boundary conditions, and tak- 
ing the two Ising spin configurations from the upper half 
and lower half of the cylinder. For XOR-Ising inter- 
faces obtained in this way, the left-crossing probabilities 
For the measured XOR-Ising loop density in [Table | the are consistent with SLE 4 1/v ^-i l/Va-r This is some " 

what surprising, since these boundary conditions satisfy 
the domain Markov property at the discrete level, so 



appropriate value of x is y/2; this value is also consistent 
with the variance and gasket dimension measurements. 

We generated many XOR-Ising loop configurations on 
2 7 x 2 15 cylinders, and using Marshall's Zipper program 
||36f computed the electrical resistances involving the non- 
contractible loops. As shown in Figures || and ||, there is 
excellent agreement between the measured distributions 
and the distributions for contour lines of the GFF with 
height spacing y/2 times as large as for CLE4. 



XOR-ISING INTERFACES IN 
DOMAINS WITH BOUNDARY 

To further test the connection between XOR-Ising in- 
terfaces and the Gaussian free field, we measured left- 
crossing probabilities of the XOR-Ising paths connecting 





two boundary points of a domain with boundary. Fig- 
ure [| illustrates the coordinate 9 against which we plot- 



ted left-crossing probabilities in Figure 5 There are a 



variety of boundary conditions for the two Ising mod- 
els that force an interface in the XOR-Ising configura- 



tion. For the solid boundary conditions in Figure 1, the 



left-crossing probabilities of the interface are consistent 



FIG. 4. (Color online) O n the left is the XOR-Ising interface 
from bottom to top from Figure 1; on the right is the corre- 
sponding XOR-Ising interface from left to right. When the 
domain is conformally mapped to the upper half plane with 
start- and end-points of the interface mapped to and 00, the 
green curves get mapped to rays re 1 for con stant 9. When 
plotting left-crossing probabilities in Figure gj , for the region 
on the left we took sites on the short diagonal and compared 
the site's 8 with the site's left-crossing probability; for the 
region on the right we did this for sites on the long diagonal. 
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Pr[left] Three left-crossing probability curves: 
1 1) XOR-Ising interface with 
black-white / all-black b.c. 

a) connecting near opposite 
corners of rhombus 

b) connecting far opposite 
, corners of rhombus 



2) SLE 4 y2_i,V2-i 




Pr[left] Three left-crossing probability curves: 
1 1) XOR-Ising interface with 

alternating / alternating (with phase shifts) b.c. 

a) connecting near opposite corners of rhombus 

b) connecting far opposite corners of rhombus, 
2) SLE 4 ,- 




Pr[left] 

1 r 



Two left-crossing probability curves: 

1) XOR-Ising interface from cylinder 
with antiperiodic b.c. 

2) SLE 4a/ ^ / 2_ 11/ ^2_ 




FIG. 5. (Color online) Left-crossing probabilities for XOR- 
Ising interfaces with different boundary conditions and differ- 
ent geometries, compared with the left-crossing probabilities 
for SLE4 iP ,p for different p's. (The top, middle, and bottom 
plots show 3 curves, 3 curves, and 2 curves, respectively; in 
each plot the curves are essentially indistinguishable.) The ge- 
ometry of the domain appears not to matter, consistent with 
conformal invariance of the XOR-Ising interfaces. Different 
boundary conditions correspond to different p's. The left- 
crossing probabilities for SLE4 iPiP were measured by sampling 
the left-crossing probabilities of a contour line in a Gaussian 
free field with ±(1 + p)\ boundary conditions p]j |. 



one might expect the interface to be simply SLE4. But 
evidently these boundary conditions act differently on 
smooth boundaries than they do on rough boundaries. 

We thank Scott Sheffield, Gady Kozma, and Bernard 
Nienhuis for useful discussions. 
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